CLAIMS 

I claim: 

1. A method for processing data to provide a forewarning 
of a critical event, comprising: 

acquiring a plurality of sets of data for at least one 
channel of data for at least one test subject or process; 

computing a renormalized measure of dissimilarity from 
distribution functions derived from a phase space for each 
respective channel of data; 

comparing said renormalized measure of dissimilarity to a 
threshold {Uc) for a number of occurrences {Nqcc) to indicate a 
condition change in said renormalized measure of 
dissimilarity; 

detecting a simultaneous condition change in a plurality 
(Nsim) of renormalized measures of dissimilarity to determine a 
forewarning of the critical event; and 

wherein said one channel of data corresponds to a parameter 
that is calculated from a plurality of parameters 
corresponding to a plurality of channels of data. s 

2. The method of claim 1, wherein the test subject is a 
human patient . 

3. The method of claim 1, wherein the test subject is a 
mechanical device or physical process, 

4. The method of claim 1, wherein the process-indicative 
data is three-phase electrical power. 

5. The method of claim 1, wherein the process - 
indicative data is vibration mechanical power. 

6. The method of claim 1, wherein the process- indicative 
data is a difference between two channels of EEG data. 
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7. The method of claim 1, further comprising: 
performing a first filtering of each set of data with a 

zero-phase quadratic filter that filters out high-frequency 
artifacts; and 

performing a second filtering of each set of data with a 
zero-phase quadratic filter to filter out low-frequency 
artifacts . 

8. The method of claim 1, further comprising: 

sorting the data values into ascending order from a 
minimum to a maximum; 

determining the number of unique signal values in) and 

the corresponding relative occurrence frequency (Fj^) for each 

unique signal value (V;,) ; 

displaying a graph of frequency {Fj^) versus values (v^) in 
each bin in a connected phase space; and 

discarding data that has [Vj^ - (N/n)]/a^ > z, where the 

value of z is determined by solving l/n = M erfc(2/V2), and 
is the standard deviation in the occurrence frequency. 

9. The method of claim 1, with an alternative embodiment 
for event forewarning, comprising determining a sequence of 
renormalized phase space dissimilarity measures from data sets 
for the test subject or process; summing said renormalized 
measures into a composite measure, C^, for the i-th data set; 

performing a least -squares analysis over a window of m points 

of the said composite measure to obtain a straight line, 
y^=ai+jb, that best fits said composite data in a least-squares 

sense; determining the variance, Gi^ = Si {y± - C±)^/(m-l), of 
said composite measure with respect to the straight line fit; 
obtaining the variability of the sequel window of m sequential 

points via the statistic, G = Si (yi - Ci)^/ai^; comparing said 
value of G to the running maximum value of the same statistic, 
Gmax/ determining the forewarning of or failure onset of a 
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critical event (such as a machine failure) , when G is 
significantly more than Gmax; obtaining the ratio, R - 
(^max) k/ ^^maJ k-if ^f the present and previous running maximum in 
G; and determining the forewarning of a critical event when R 
is significantly more than some limit, 

10. A method for processing data to provide a 
forewarning of a critical event, comprising: 

acquiring a plurality of sets of data for at least two 
channels of data for at least one test subject or process; 

computing to a multi-channel time-delay phase-space (PS) 

construction, which has the form: y(i) = [s(l)^, 3(1) # 

s(l)^^2^,... , s(2), s(2)^^x, s(2)^^2^, s(c)i, s(c)^^x, s (c) ^^2^, ..J], 

where s(c) denotes the symbolized data for c-th channels- 
computing a renormalized measure of dissimilarity from 

distribution functions derived from the (non) connected phase 

space for the multi -channel of data; 

comparing said renormalized measure of dissimilarity to a 

threshold {Uc) for a number of occurrences (Nqcc) to indicate a 

condition change in said renormalized measure of 

dissimilarity; and 

detecting a simultaneous condition change in a plurality 

(Nsi„) of renormalized measures of dissimilarity to determine a 

forewarning of the critical event. 

11. The method of claim 10, further comprising: 
performing a first filtering of each set of data with a 

zero-phase quadratic filter that filters out high-frequency 
artifacts; and 

performing a second filtering of each set of data with a 
zero-phase quadratic filter to filter out low- frequency 
artifacts . 

12. The method of claim 10, using an alternative 
embodiment for event forewarning, comprising determining a 
sequence of renormalized phase space dissimilarity measures 
from data sets collected during increasingly severe fault 
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conditions; summing said renormalized measures into a 
composite measure, for the i-th data set; performing a 

least -squares analysis over a window of w points of the said 
composite measure to obtain a straight line, yi = ai+jb, that best 
fits said composite data in a least-squares sense; determining 
the variance, Oi^ = Si (yi ~ Ci)^/(7n-l), of said composite 
measure with respect to the straight line fit; obtaining the 
variability of a sequel window of m sequential points via the 

statistic, G = Si {y± - Ci)^/ai^; comparing said value of G to 
the running maximum value of the same statistic, Gmax; and 
determining the onset of a critical event, such as forewarning 
of a machine failure, when G is significantly more than G(non- 
end-of -life) , or when R is significantly more than i?(non-end- 

of-life) , or detection of failure onset when G is 
significantly greater than G(end-of -lif e) . 

13. The method of claim 10, wherein the test subject is 
a human patient. 

14. The method of claim 10, wherein the test subject is 
a mechanical device or physical process. 

15. The method of claim 10, further comprising: 
sorting the data values into ascending order from a 

minimum to a maximum; 

determining the number of unique signal values (n) and 
the corresponding relative occurrence frequency (Fj^) for each 
unique signal value {Vj^) ; 

displaying a graph of frequency (F^) versus values (Vj^) in 
each bin in a connected phase space; and 

discarding data that has [v^, - {N/n)]/(j^ > z, where the 

value of z is determined by solving 1/n = H erfc(z/V2), and 
is the standard deviation in the occurrence frequency. 
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16. The method of claim 10, wherein the process- 
indicative data is three-phase electrical power. 

17. The method of claim 1, wherein the process- 
indicative data is vibration mechanical power. 

18. The method of claim 1, wherein the process- 
indicative data is a difference between two channels of EEG 
data. 
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